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(54) Title: METHOD FOR DETERMINING HAEMODYNAMIC INDICES BY USE OF TOMOGRAPHIC DATA 
(57) Abstract 

Haemodynamlc indices of an organ or a part of 
tissue are determined from a time scries of tomographic 
data obtained by means of Magnetic Resonance Imaging. 
Maps of indices are produced, being significant of the 
dynamics of the capillary tissue flow acquired during rapid 
bolus injection of a tracer that stays mainly intravascular. 
The method may be used for evaluating the efficacy of 
a drug on an organ, or for obtaining information of the 
likelihood of recovery of an organ or part of tissue upon or 
during a period of insufficient vascular supply or during the a 
progression of a chronic disease. The method may be used s 
for discriminating between relevant therapy of an organ. 
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Instead, image elements containing partly tissue, partly blood vessels, must be usied in 
order to characterise the mere shape of the arterial input curve to the tissue. In such 
cases, absolute values of flow and volume cannot be found, and therefore a normalisation 
routine allowing (a) comparison of serial measurements in a single subject, (b) comparison 
5 among subjects and (c) absolute quantification the case of susceptibility MRI of the brain is 
described below. 

OVERVIEW 

10 

The techniques consist of the following steps: 

A. Conversion of tomographic images into images representing concentrations of contrast 
agent as a function of time. 

B. Identification of relative or absolute (i.e. in units identical to those of the tomographic 
15 images after the above mentioned conversion) arterial tracer concentrations from the 

image data. 

C. In the case of relative arterial tracer concentrations. 

normalisation of the arterial input area to the injected dose of contrast agent per kg 
body weight. 

20 D. Optional correction of tissue tracer curves for delays. 
E. Determination of 

(i) absolute or relative tissue blood flow 

(ii) tissue impulse response function 

by deconvolution of tissue concentration time curves by the arterial tracer concentration 
25 in each image element. 

P Determination of tissue blood volume by determining the area under the tissue first- 
pass concentration curve. 

G. Determination of tissue mean transit time by the tissue blood volume - blood flow ratio. 

H. In the case of MR susceptibility contrast imaging using Gd-chelates in brain tissue. 
30 conversion to absolute blood flow and blood volume by a pre-determined constant. 

I. Determination of the distribution of flow or transit times in each image element from the 
residue function (nomialised impulse response function) determined In E. 

J. Comparison of distributions of relative flows to a predetermined distribution found for a 
normal organ; here brain. 
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K. Quantification of the distribution of flows in terms of the extraction fraction of a given 
solute with specified capillary permeability, in cases where imaging is performed with 
MR imaging with microvascular weighting, or microvascular volume can otherwise be 
inferred from total blood volumes. 

5 

1 . Conversion of tomographic data into tracer concentration images 

The method and associated software handles various modalities, depending on whether 
tracer injection changes signal intensity from baseline in a linear or logarithmic fashion 
10 upon tracer arrival. With a specified option, acquired tomographic images during tracer 
bolus injection are converted into concentrations as a function of time, with measurement 
time points spaced equally in time. 

1.1. In the case of MR images, weighted towards the transverse relaxation times (Tz or T2'), 
15 (typically acquired In brain), the tissue concentration as a function of time, Ci(t), is typically 
obtained by the fonnula 



30 



C,(0 = -/flog 



f S{t) 

[S{to)) 



ITE Eq. 1 



20 where S(to) is the signal intensity before contrast injection (formed as the average of signal 
intensities up to tracer arrival), S(t) is the signal intensity at time t, and TE is the echo time 
used in the sequence. Here, /c is a constant characteristic of the tissue and contrast agent. 

1.2. In the case of MR images weighted towards the longitudinal relaxation time (T,) 
25 images (typically acquired for intravascular tracers in the heart), concentrations are 

assessed directly by the change in longitudinal relaxation rate, ARi, as derived from signal 
intensity changes 



C,(0 = W A/?i Eq. 2 

where 9? is the relaxivity of the applied contrast agent. 



1.3. In the case of Computed Tomography (CT) images, concentrations are assessed by 
the change in image intensity, measured in Hounsfield units, AH 
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C^(0 = a: AH Eq. 3 

where k is the characteristic X-ray absorption of the contrast agent. 

5 

2. Identification of Arterial Vessels 

For identification of arterial vessels in the image, the algorithm provides the choice of first 
producing an image that guides this process (2.a.), or proceeding directly to calculation of 
10 tissue flow (2.b.) 

2.a. If so specified by an option, the concentration time curve of each image element is 
fitted to a gamma variate function by nonlinear least-squared regression (by means of a 
fletcher-algorithm) over the time interval until visible tracer re-circulation occurs (the re- 
15 circulation is specified by the user, determining 

(a) Tracer arrival time 

(b) Area under first pass (this quantity is proportional to the local blood 
volume). 

These two quantities are stored as floating point binary image file, allowing visualisation as 
20 an image where each image pixel corresponds to quantities (a) and (b), respectively. An 
artery is then visually identified in these images by 

(a) Anatomical location. 

(b) Early tracer arrival relative to the arrival in tissue. 

(c) Large area under first pass (corresponding to a large blood volume and 
25 thereby a large portion of the vascular volume being contained in a voxel). 



2. b. By observing images of tracer concentration as a function of time in a tool that allows 
visualisation of the time-course of single pixels, arterial levels can be identified. 

30 Upon identification of the arterial input function, the con^esponding signal intensity time 
curve is fed to the program below in the form of an ASCII file. 

3. Normalisation of Arterial Input to Allow Comparison of Serial Measurements 



